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Abstract 

Small chemical sensors are subjected to adsorption-desorption fluctuations which 
usually considered as noise contaminating useful signal. Based on temporal proper- 
ties of this noise, it is shown that it can be made useful if proper processed. Namely, 
the signal, which characterizes the total amount of adsorbed analyte, should be sub- 
jected to a kind of amplitude discrimination (or level crossing discrimination) with 
certain threshold. When the amount is equal or above the threshold, the result of 
discrimination is standard dc signal, otherwise it is zero. Analytes are applied at low 
concentration: the mean adsorbed amount is below the threshold. The threshold is 
achieved from time to time thanking to the fluctuations. The signal after discrimi- 
nation is averaged over a time window and used as the output of the whole device. 
Selectivity of this device is compared with that of its primary adsorbing sites, based 
on explicit description of the threshold-crossing statistics. It is concluded that the 
whole sensor may have much better selectivity than do its individual adsorbing sites. 
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1 Introduction 

Detectors of chemical substances are usually based on selective adsorption- 
desorption (binding-releasing) of analyzed chemicals by specific adsorbing sites 
(receptor molecules). The receptor molecules are attached to an electronic 
device able to measure the amount of the analyte adsorbed during the binding- 
releasing process. The device may be either a MEMS device, such as quartz 
crystal microbalance [1,2], or vibrating/bending cantilever [3], or field effect 
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transistor [4], or other [5]. The device with the receptor molecules is called 
chemical sensor or detector. In order to be useful, the detector must be able 
to discriminate between different chemicals, to be selective. Its selectivity is 
normally the same as that of its receptor molecules (see Eqs.(6,7)). 

The size of industrial sensors has constant tendency to decrease [3]. The power 
of useful signal produced by a small detector becomes very small. As a result, 
noise of the detector itself constitutes a substantial portion of its output signal. 
Depending on its construction, there are several reasons for a small detector 
to be noisy [6]. One type of noise is due to the fact that the adsorption- 
desorption process is driven by brownian motion, which is stochastic. As a 
result, the instantaneous total amount of adsorbed analyte is subjected to 
irregular fluctuations visible in the output signal. This noise is called the 
adsorption-desorption noise [7]. It is present in any small detector which is 
based on binding-releasing of analyte. The adsorption-desorption noise can 
dominate over all other types of intrinsic noise [8]. 

In this paper only the adsorption-desorption noise is taken into account. The 
detector is expected to be a threshold detector (ThD), Fig.l. 



>- 



Soit) 




m 




S 




ThU 


11 m nil II 


TAU 
















- PSU 











A R 



Fig. 1. Schematic picture of threshold detector. A — analyte molecules; R — adsorp- 
tion sites; PSU — primary sensing unit; ThU — threshold unit; TAU — temporal 
averaging unit. 

Namely, the fluctuating signal characterizing the amount of adsorbed analyte 
in the primary sensing unit (PSU in Fig.l) is fed into amplitude discriminator 
unit (threshold unit, ThU in Fig.l). The threshold unit is characterized by a 
certain threshold. It has zero as its output if the adsorbed amount is below 
the threshold, and it outputs standard dc signal while the adsorbed amount 
is equal or above the threshold. The output of ThU is averaged over a sliding 
time- window to have final output practically time- independent. This signal is 
considered as the output of the ThD. 

In this paper, the temporal properties of the binding-releasing stochastic pro- 
cess are utilized to characterize the outputs of ThD if two analytes are sep- 
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arately presented at equal concentrations. This allows to compare selectivity 
of TliD with that of its receptor molecules. The main conclusion is that the 
ThD may be much more selective than do its adsorbing sites. 



2 Definitions and assumptions 

The adsorption-desorption process is described by the following association- 
dissociation chemical reaction 

A + R^ AR, (1) 



where A, R, AR denote molecules of analyte, adsorption site or receptor, and 
analyte-receptor binary complex, respectively. At constant temperature, the 
rate constants, k^, /c_ are time-independent. They can be determined either 
from experimental measurements, or estimated theoretically [8]. Let N denotes 
the total number of receptor molecules per detector. The analyte is presented 
at concentration c. The probability p for any R to be bound with A isQ 



k+c 

P 



k+c + k. 



The mean number of adsorbed molecules, (ra), can be calculated as follows: 

(n) = pN. 



If two different analytes Ai, A2 are tested at the same concentration, either 
the Eq.(2), or experimental measurements will give two values, pi, P2- We say 
that the receptor molecule has selectivity with respect to Ai, A2, if pi 7^ P2 



see [9], where Eq.(2) is justified. 
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(expect, pi > J92). The molecular selectivity, fi, is defined asQ[3] 




(3) 



The primary signal, So{t) in Fig.l, usually increases if the number n of ad- 
sorbed molecules increases: 



where the exact dependence of So on n is determined by the sensor construc- 
tion and the transduction mechanism it employes. For simplicity, it is expected 
that in the case of gravimetric sensor, Ai and A2 have equal molecular masses. 

Define selectivity 6 for a whole detector in terms of final output signal {S in 
Fig.l) as follows: 



where 5*1, 5*2 are the final outputs for analytes Ai, A2, respectively. 

Both 5*0 (t) and n{t) are subjected to adsorption-desorption noise. In a detec- 
tor without the threshold unit, the final output signal can be made linearly 
proportional to the mean number of adsorbed molecules: 



This is achieved either by temporal averaging, or by choosing large detector 
with powerful primary signal in which contribution of adsorption-desorption 
fluctuations is not visible. Substituting (6) into (5) one obtains for selectivity 



^ If one do not expect that pi > p2 than Eq.(3) should be replaced by /u = 
\ln{pi/p2)\. 

^ This definition of selectivity differs from used in chemistry the specificity of as- 
sociation which is expressed in terms of dissociation constant. For analyte A, the 
dissociation constant is defined as [A]i/2 = k^/k+. Eq.(2) can be rewritten using the 
dissociation constant: p = 1/(1 + [^]x/2/c). Prom this equation it is clear that an- 
alytes with different dissociation constants have different binding probabilities and 
vice versa. This proves suitability of both descriptions, even if numerical values of se- 
lectivity expressed in terms of dissociation constants, say as = ln([742]i/2/[^i]i/2)) 
will differ from used here. The fi values can be expressed in terms of dissociation 
constants: fi = ln((c-|- [yl2]i/2)/(c + [^i]i/2))- The main difference between the fj,' 
and )U is that the latter depends on concentration. This is in accordance with situ- 
ation in natural olfactory systems where discriminating ability usually depends on 
concentration [10]. 



n > n' ^ Sq > S'q 



(4) 




(5) 



Si ~ PiN, 1 = 1,2. 



(6) 
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of a conventional detector 

5 = In — - = II. 
P2N 



(7) 



Thus, selectivity of detector in which the fluctuations are averaged out either 
immediately after the primary sensing unit, or inside it is equal to that of its 
individual adsorbing sites. 

The threshold unit, ThU, rises a threshold which the So must overcome in 
order to make possible further stages of processing. The crossing may happen 
from time to time thanking to the adsorption-desorption fluctuations. Due to 
(4) , the threshold can be characterized by the number A^o of analyte molecules 
which must be adsorbed before the nonzero signal appears at the output end 
of the ThU. It is assumed that the ThU is ideal in a sense that the Nq is the 
exact value which is not subjected to fluctuations. If A'o is achieved, the ThU 
has standard constant signal as its output. The signal does not depend on the 
exact value of n{t) provided it is above or equal to A^o- 

Denote by T the temporal window over which the averaging is made in the 
TAU (Fig.l), and by Tb, Ta (Tb + Ta = T) the total amount of time during 
which n{t) is below or above the threshold, respectively, when < t < T. The 
flnal output, S in Fig.l, should be linearly proportional to T^/T. This gives 
for the selectivity of ThD: 

<^ = ln^, (8) 
where Tqi, Ta2 correspond to Ai, A21 respectively. 



3 Estimation of selectivity 



In accordance with (8), it is necessary to estimate the total amount of time the 
n{t) spends above the threshold when t e [O;!"]. This can be done by adding 

together lengths of all separate intervals during which n{t) > Nq continuously. 
Denote by ( the number of those intervals, and by T^, 1 < k < the length 
of the k-th continuous interval. Then 

rp _ \ " rpk _ A_ \ " rpk _ rp 1<^'<C 

i<fc<C ^ i<fc<C C.Ir'^y^ o-' 

l<fe<4 

where is the length of k-ih. continuous interval during which n{t) < Nq. If 
T together with c, /c_ ensures that C, is large, then the last expression can 
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be rewritten in the following form 

T„ = T , (9) 

be ' ac 

where The, Tac are the mean lengths of the continuous intervals. For the The, Tac 
the following expressions have been obtained [11] based on the Kolmogoroff 
(or backward Master) equation: 

If two analytes, Ai, are considered, then in (10), (11), fc- and p should be 
replaced with pj, i = 1,2, respectively. Substituting (10), (11) into Eq. 
(9) one obtain^ 

T,=T y: cya-pf-'. (12) 

No<l<N 

Considering (12) for two analytes, use it in Eq.(8). This gives 

r. , N0<l<N . X 

No<l<N 

The last equation can be replaced by a transparent estimate if one use the 
following inequality 

> f ^ ] " (14) 

No<l<N 

which is proven in [11]. Substitution of (14) into Eq.(13) gives 

6 > iV ^° where po = No/N. (15) 

l-pi 

Taking into account that the total number of adsorbing sites, A^, as well as A^^o 
can be very large, it is clear from the estimate (15) that 6 can be much larger 



The following relation is used: E C'jyp (1 — p) = 1. See also [9], where 

0<1<N 

equivalent to (12) conclusion is obtained based on simpler and less rigorous reasoning 
not using Eqs. (10), (11). 
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than /i, provided the fraction (po — — Pi) is not very smalL It must be 

at least positive, which requires 

po > pi, or piN < No. (16) 

Taking into account that pi increases with concentration (sec Eq.(2)), inequal- 
ity (16) can be considered as imposing an upper limit for concentration c at 
which the effect of selectivity improvement might be expected based on the 
estimate (15). It is worth to notice that when condition (16) holds, the mean 
amount of adsorbed analyte is below the threshold one, and threshold crossing 
may happen only due to fluctuations. 



4 Numerical examples 

As one can conclude from the estimate (15), the selectivity improvement is 
higher for higher A^"o. On the other hand, one cannot chose the Nq as high 
as desired because the ThU in Fig.l is expected to be ideal. If one chose 
Nq = 100 then the ideality means that the threshold level in the ThU is 
allowed to have less then 1% jitter. Similarly, if one chose A^o = 10^ then 
the threshold level must be kept with better than 0.01% precision. Otherwise, 
noise in the threshold level should be taken into account in the reasoning of 
n.3, and this will lead to a less promising estimate. Another conclusion, based 
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Table 1 

The rate constants used in the examples of Table 2 and in Fig. 2. 

on the estimate (15), suggests that the smaller is the concentration (smaller pi) 
of the analytes, the better is discrimination between them. But in this case the 
threshold will be achieved during small fraction of time spent for measuring. 
As a result, the output signal will be very small and may be lost in the TAU 
unit. It is natural to require that the output signal for more affine analyte 
is higher than the 10% of the maximal output signal, which is produced if 
n{t) > Nq all the time. Taking into account Eq.(12) this leads to the following 
constraint 

ri= E C^Wi(l-Pi)^~'>0.1. (17) 

One more constraint comes from assumption of large ( which is made for 
derivation of Eq.(9). If the measuring (averaging) time T is to be short enough. 
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0.05 


18 


0.16 


4 • 10-5 



Table 2 

Numerical examples of improved selectivity. The rate constants for the analytes are 
shown in the Table 1. (5 is calculated here by means of the exact expression (13), ri 
— as shown in Eq.(17). 




9e-05 9.6e-05 0.000108 9e-05 9.6e-05 0.000108 



Fig. 2. Concentration dependences of selectivity for the examples of Table 2. Con- 
centration (x-axis) is given in M. The A'^, A^o values in a and b correspond to the 
first and second rows of Table 2, respectively. The Sest corresponds to the right hand 
side of the inequality (15). 

say T = 1 s, then the mean frequency of crossing the threshold should be 
high enough in order to have, e.g., ( > 1000. This could be achieved if the 
mean durations of being continuously above and below the threshold are short 
enough. If fabc = + Tac, then Eqs.(10),(ll) give 

The ( > 1000 could be ensured by the following inequality 

T/f^bc > 1000, 

or, choosing T = 1 s, 

fabc < 0.001 s. 

Two examples satisfying this constraints are shown in the Table 2. Concen- 
tration dependencies of n, 5, and the estimate (15) are shown in Fig. 2. A short 
segment of the trajectory n{t) modelled on PC is shown in Fig.3. 

5 Conclusions and discussion 

In this paper, selectivity of chemical sensor is compared with that of its pri- 
mary receptors (adsorbing sites). The sensor is expected to be a small one, 
in which the main source of noise is due to the adsorption-desorption fluc- 
tuations. In the sensor considered, the signal from the primary sensing unit 
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Fig. 3. Short segment of the trajectory n{t) modelled on PC for the Example 1 of 
Table 2. Time (x-axis) is given in seconds. 

is immediately subjected to the amplitude discrimination defined in the In- 
troduction, and obtained piecewise-constant signal {L{t) in Fig.l) is averaged 
over a time window. The averaged signal {S in Fig.l) is taken as the output 
of whole sensor. 

The threshold-crossing statistics derived from the exact description of the 
adsorption-desorption stochastic process is used for estimating selectivity. As 
a result, it is concluded that selectivity of this sensor can be much better than 
that of its primary receptors. The effect may be expected in a limited range 
of concentrations of analytes, which depends on the threshold level. For high 
concentrations the selectivity falls to that of the primary receptors (Fig. 2), and 
for low ones the output signal will be too small even for more affine analyte. 
The best situation is expected when the mean number of bound receptors 
is just below the threshold one, and the threshold is frequently crossed due 
to the presence of fluctuations. Thus, in practical realization a possibility of 
tuneable threshold should be considered. 

Usually, noise in sensory devices is taken as unfavorable factor[^. In this con- 
sideration, the presence of noise looks like factor improving the sensor per- 
formance. But with the ideal threshold unit in hands much can be done even 
without noise. Expect that the noise is initially averaged out either by spatial 
averaging (choosing big primary unit with large A^), or by temporal averaging 
(interchanging TAU with ThU in Fig.l). The averaged signals for the Ai, A2 
can be very close (see Eq. 6), but the ideal ThU with tunable threshold will 
be able to discriminate perfectly between them. Thus, even if the fluctuations 
in this sensor are made working, the answer what is better to do first for 
the practical purposes: the amplitude discrimination, or temporal averaging, 
depends on physical parameters of the environment in which the sensor oper- 
ates, and on physical characteristics of the sensor itself, including intensity of 
noises other than the adsorption-desorption one. Interesting, in natural olfac- 
tory systems, a kind of amplitude discrimination is made immediately after 
the primary reception [11,13]. Also in those systems the threshold is tunable 
due to adaptation of individual neurons. 



^ but see [12], where some characteristics of noise are employed for discriminating 
purposes. 
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